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Abstract 

We present a detailed analysis of resonance contributions in the context of higher 
twist effects in the moments of the proton spin structure function g\. For each of 
these moments, it is found that there exists a characteristic Q 2 region in which 
(perturbative) higher twist corrections coexist with (non-perturbative) resonance 
contribution of comparable magnitude. 
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1 Introduction 



High energy lepton scattering is a well established tool to investigate the structure of 
the nucleon. We restrict ourselves to charged leptons (e or /i); the exchanged virtual 
photon transfers four- momentum = (qo,q), with the resolution determined by the 
virtuality Q 2 = —q 2 = q 2 — q 2 . At Q 2 ^> lGeV 2 deep-inelastic scattering resolves the 
partonic constituents (quarks and gluons) of the nucleon. At Q 2 < lGeV 2 , on the other 
hand, the excitation of nucleon resonances and multi-pion continuum states is important. 
Exploring the transition between partonic and hadronic scales is of great significance to 
our understanding of the nucleon. The aim of the present paper is to discuss polarized 
lepton-nucleon scattering in kinematic regions where both hadron and parton degrees of 
freedom are expected to coexist. 

The response of the nucleon is expressed in terms of the hadronic tensor 

W„ u (x,Q 2 ) = ^J2( 2 ^( p +Q- p x)(N(P,S)\U0)\X(Px,Xx)} 

47T x 

(X(P X ,\ X )\J V {Q)\N(P,S)) (1) 
= W {S) + W {A) 

The matrix elements of the electromagnetic current describe the transition of a nu- 
cleon with four-momentum P, invariant mass M (P 2 = M 2 ) and spin S to a hadronic 
final state X with four- momentum P x and polarization Ax- The sum in (|I|) implies an 
integration over three-momentum, ^ 2 w)^2P X o ' anc ^ ^ e normauza ti° n of \N) and \X) is 
(N(P', S)\N(P, S)) = 2P (2n) 3 5 3 (P' - P)5 S>S , 

The symmetric part Wffl involves the spin independent structure functions F 12 measured 
in the scattering of unpolarized particles. The antisymmetric term, 



W$ = ie, vXa q x 



^Q%, + 9^Q 2 ) (P . qS ,_ q . sn 



(2) 



P-q (P-q) 2 

introduces the spin structure functions g\ and g%. The nucleon spin vector is 5°" = 
\u(P, 5 , )7 <t 7 5 m(P, S) with Dirac spinors normalized as uu = 2M. The structure functions 
depend on the Bjorken variable x = Q 2 / (2P ■ q) and on Q 2 . 

Spin structure function data have been taken at SLAC, CERN and DESY |1|, 0, |3|, f|, 
primarily in the partonic high Q 2 range. Polarized deep-inelastic scattering in the reso- 
nance region was measured by the E143 collaboration at SLAC In the first part of 
our study we combine these data with other available information from the photo- and 
leptoproduction of nucleon resonances and investigate their contribution to the moments 
of the proton spin structure function g^. 

The influence of resonances and non-resonant low-mass excitations turns out to be quite 
significant for Q 2 < 4 GeV 2 , as we shall demonstrate. For example, at Q 2 = 2 GeV 2 
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they account for as much as 20% of the first moment of g±. Similar observations have 
been made for unpolarized deep-inelastic scattering ||. 

In the second part we use the QCD operator product expansion and extract twist-4 matrix 
elements from the leading moments of g\. For the first moment such an analysis has been 
carried out in great detail in ref . . We find substantial higher-twist contributions to the 
first, third and fifth moments of g\ for Q 2 < 2,4 and 10 GeV 2 , respectively. We examine 
target mass effects and investigate the different components of the higher-twist pieces of 
gi . It turns out that contributions from elastic scattering, low-mass hadronic excitations 
and the partonic high-mass continuum are all of similar importance. We comment on 
the applicability of the twist expansion and recall basic ideas of parton-hadron duality. 
Altogether our results emphasize the need for high-precision experiments in the resonance 
region, to be performed at the Jefferson laboratory 0. 



2 Twist expansion of g 1 

In this section we briefly summarize results from the operator product expansion for the 
nucleon spin structure function g\ (for details see e.g. [§). Following the conventions of 
ref. || we introduce the n-th moment of g\ as: 

gP(Q 2 ) = [ tixx^gAx^Q 2 ) (with 11 = 1,3,5...)- (3) 
Jo 

Note that the upper limit of integration includes the contribution from elastic scattering. 
Its presence results from the fact that the operator product expansion, applied to deep- 
inelastic scattering, implicitly involves a sum over all final hadronic states including the 
nucleon itself. The importance of the elastic component in a QCD analysis of structure 
function moments has been emphasized in ref.0. 

At large momentum transfers, Q 2 ^> Aq CD , the moments (Q) can be written in terms of 
the twist expansion [jiCfl : 

g{ n \Q 2 )= E ^M 1 - (4) 



=2,4,. 



"Twist" is a useful bookkeeping device to classify the light cone singularity of the coeffi- 
cients in the QCD operator product expansion. Let a local operator in this expansion be 
a Lorentz tensor of rank r with (mass) dimension d, and let a < r be the "spin" associated 
with this operator. Then twist is defined as r = d — a. The functions /^(Q 2 ) are re- 
lated to nucleon matrix elements of quark and gluon operators with maximal twist r. 
Their leading (logarithmic) Q 2 -dependence can be calculated perturbatively as a series 
expansion in the strong coupling constant a s . It should be mentioned that, due to the 
asymptotic nature of QCD perturbation theory, a systematic separation of the twist ex- 
pansion and the perturbation series for is non-trivial and still a matter of ongoing 
investigations (for detailed discussions see e.g ref.pl]]). 
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Up to corrections of order 1/Q A one finds 



9 { r\Q 2 ) = ^.-i(Q 2 ) + ^^^(^an + i(Q 2 )+4 4 + i(g 2 )) (5) 

AM 2 




The coefficients a n represent the genuine twist-2 contributions to g[ n \ They depend only 
logarithmically on Q 2 and dominate for Q 2 much larger than a typical hadronic scale, say 
the squared nucleon mass M 2 . The second term in (|5]) arises from target mass corrections 
| 12| . They are determined by the twist-2 pieces a n and the twist-3 corrections d n related 



to moments of the spin structure function g 2 [13] : 

dn-i = 2^ n) + ^-/yV" - 0[ — \. ((>) 




The true twist-4 contributions in eq.(^J) are denoted by f n +i. For higher moments, n > 1, 
several matrix elements of twist-4 are involved. Their sum gives the coefficient f n+ \ in 



In our work twist-2 contributions are defined through moments of presently available NLO 
parametrizations of g\ |14| . The extraction of higher twist contributions from structure 



function data has been a subject of recent studies 0. The active interest in these quan- 
tities derives from the fact that they are related to matrix elements which are sensitive 
to quark-gluon interactions in the nucleon. For example, one has || |T5|| : 

2f 2 (Q 2 )M 2 S» = J2e}(N(P,S)\gj> f G^^ f \N(P,S)). (7) 
/ 

The sum is taken over all quark fields ipf with flavor / and charge e/, and G^ v stands for 
the dual gluon field strength tensor (g denotes the QCD coupling strength). 



3 Helicity amplitudes 



In this paper we investigate contributions to the proton spin structure function gf resulting 
from the electro-production of nucleon resonances, as well as from the production of 
continuum states in the deep-inelastic regime. Resonance contributions are conveniently 
described in terms of helicity amplitudes |Tj| : 



G r . 



2M 



X{P X ,X' 



m 



J(0) 



N(P,X 



(8) 



We choose q/\q\ as the spin quantization axis. The amplitude G m represents the produc- 
tion of a hadronic state X with spin projection A' following the absorption of a virtual 
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photon with polarization (helicity) m = ±1, on a nucleon with spin projection A = —1/2. 
The photon polarization vectors are e 11 = (0, +1, —i, 0)/y/2 and e° = (|q|, 0, 0, u)/Q, with 



Combining eqs.(|T]) and (§) gives: 
1 



9i 



1 + 



M 2 5(W 2 - M 



la 



(9) 



with v = P ■ q/M. The final state X with invariant mass Mx has A' = +1/2 for G+, 
A' = —3/2 for and A' = —1/2 for G$. It is common to use the amplitudes (e is the 
electric charge with e 2 /An = 1/137): 



A 



M 



1/2 



W 2 -M 2 



G + , A 



M 



3/2 



W 2 — M 2 



G-, Si /2 = e 



M 



W 2 -M 2 Q 



Go, (10) 



where q* denotes the three-momentum transfer as measured in the photon nucleon center- 
of-mass frame, i.e. q* 2 = Q 2 + (W 2 — M 2 — Q 2 ) 2 /AW 2 with the total c. m. energy W. 



4 Model 

In the following we present a parametrization of the proton structure function g\ which 
is applicable at small and moderate values of Q 2 . We follow hereby closely an analysis 
of recent data in ref.[||. At small center-of-mass energies, W < 1.7 GeV, we account for 
the contribution of dominant nucleon resonances. In addition, a phenomenological non- 
resonant background is added. For large W > 1.7 GeV we use an existing parametrization 
of available data. 

The contribution of an isolated nucleon resonance to g\ is usually expressed through he- 
licity dependent virtual photon-nucleon cross sections. In terms of the helicity amplitudes 
(PH) these are defined as: 



t MTr 

^1/2,3/2 " Mr[{w _ Mr)2 ~ r y 4] 1^1/2,3/2 



a 



mt r q 2 2 

l ' 2 M R [(W-M R ) 2 + T 2 R /4}q* 21 1/2 ' ' 

LT _ MT R Q_ q* A /-, -, X 

1/2 " V2M R [(W - M R ) 2 + r|/4] W ' ' 



Here M R is the mass and T R the width of the resonance. Combining eqs. (P JlO|JTlD gives 
for the contribution of a resonance R to g\. 

2 uM-Q 2 /2 1 , Q lt\ / 19 x 

= 4 ^a 1 + QV^ 2 ( 2 + ^ 2 J' ( } 
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where the photon-nucleon cross sections refer to the excitation of R. At low W the helicity 
amplitudes are reasonably well known only for the photoproduction of the prominent 
nucleon resonances. In the case of electro-production accurate data are rare (for a review 
see e.g. |T7|]). We restrict ourselves to the dominant low mass resonances A (1232), 
Sii(1535), and -D 13 (1520). Our parametrizations of the corresponding helicity amplitudes 
are summarized in eqs.([I3],[n]), with parameters given in Table 1. 

At low center of mass energies the excitation of the A(1232) resonance is of particular 
importance. At small Q 2 it is dominated by a magnetic dipole transition which implies 



A3/2/A1/2 « v^- Indeed, for real photons one finds ^3/2/(^^1/2) ~ 1-064 |R§. At 
large momentum transfers, Q 2 ^> 1 GeV 2 , perturbative QCD gives A 3 / 2 /Ai/ 2 ~ 1/Q 2 - 
However, it has been observed that even at Q 2 = 3 GeV 2 the magnetic dipole transition 
still dominates by far |19[ . We can therefore assume A 3 / 2 /A 1 / 2 ~ const, for Q 2 <3 



GeV . The Q -dependence of Ai/ 2 and A 3 / 2 is then extracted from an analysis of the 
Q 2 -dependence of the transverse amplitude \A?\ Jl - 

The 5*11 resonance has spin 1/2, so that the helicity amplitude A 3 / 2 is absent. We constrain 
the parametrization of the amplitude Ai/ 2 by the photo- and electro-production data from 



ref.|0|. 



For the _Di 3 (1520) the amplitude Ai/ 2 is found to be very small at Q 2 = 0. Here A 3 / 2 
dominates. On the other hand data require Ai/ 2 > A 3 / 2 for Q 2 > 1 GeV 2 |T7]. 



The parametrization in eqs. (|T3|JT^) agrees with the present, albeit quite rough, empiri- 
cal information on the Q 2 -dependence of the asymmetry Ai and the individual helicity 
amplitudes fI7] : 

12 , 1 a 1 2\ 1 / 2 _ n d ^21 



and 



\An 



1^4.1/2,3/2 I 



A 1/2 \ 2 + \A 



3/2 1 



C exp[-BQ' 



'l±Ai 



Ay I, with A\ 



l A V2 


2 " 


A3/2P 


|Al/2 


2 + 


A3/2 


2 



(13) 
(14) 



with parameters given in table 1. 





C/GeV" 1/2 


B/GcV- 2 


Ai 


T/GeV 


A(1232) 


0.293 


0.6 


-0.545 


0.12 


511(1535) 


0.07 


0.17 


1.0 


0.15 


^13(1520) 


0.16 


0.83 


1 - exp[0.3 - 1.5 g 2 /GeV 2 ] 


0.12 



Table 1: Parameters for the helicity amplitudes in eqs.(|T 
which enter in eg-fllTD, are taken from ref . 1 18 



|JI4]). The resonance widths T, 



The interference term a LT involving the longitudinal and transverse photon-nucleon am- 
plitudes is fairly unknown. Nevertheless, unpolarized scattering constrains the asymmetry 
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ratio: 

Z(T l/2 



A 2 = T 2\ < v^fo Q 2 ), ( 15 ) 

°"l/2 + a 3/2 

with R = 2a^ 2 /(a]y 2 + cr^/ 2 ). In the resonance region one finds on average R = 0.06 ±0.02 



for 1 GeV 2 < Q 2 < 8 GeV 2 and W < 1.7 GeV [pi]]. Some fraction of this value is due to 
incoherent background contributions and not related to the excitation of single nucleon 
resonances. In the following we use A 2 = 0.08. As a matter of fact, at Q 2 > 1 GeV 2 , 
A 2 contributes only very little to the structure function moments to be discussed later: 
changing A2 by 100% modifies our results for by less than 2%. 

At low energies, W < 1.7 GeV, the structure function g\ receives contributions also from 
non-resonant (multi-)meson production. However, hardly any empirical information is 
available here. We use a linear interpolation in the squared photon-nucleon center of mass 
energy W 2 which connects the inelastic threshold W 2 = (M + m^) 2 with experimental 
data at W > 1.7 GeV. 

Having modeled the structure function g\ at small center-of-mass energies, we continue 
to W > 1.7 GeV, where we use a parametrization from ref. U which reproduces data in 
the deep-inelastic region. This, finally, completes our model for g\. 

In Fig.|I| we compare our model with recent g\ data from the E143 collaboration taken 
at Q 2 = 1.2 GeV 2 . Within the admittedly large experimental errors good agreement is 
found. A comparison of gf, calculated with our model, and a parametrization of its deep- 
inelastic twist-2 part is shown in Fig.Q. At Q 2 < 2 GeV 2 , significant deviations are 



apparent. The contributions of the Su and -D13 resonances are located around x ~ 0.4 
at Q 2 = 1 GeV 2 , while the excitation of the A occurs at x ~ 0.6. As Q 2 increases 
the low mass nucleon resonance excitations become less important. Furthermore, the 
contribution of nucleon resonances moves towards larger values of x, as one can see from 
the fact that the squared invariant mass of a particular nucleon excitation is fixed at 
W 2 = M 2 + Q 2 (l — x)/x. Finally, at Q 2 = 10 GeV 2 our model coincides with the leading 
twist parametrization of ref. [14[ . 



5 Analysis of moments of structure function 

In this section we discuss the first moments of the proton structure function g\ as obtained 
from the model previously described. In particular, we investigate the importance of 
contributions from elastic scattering, resonance production, target mass corrections, and 
true higher twist. 

The elastic contribution, corresponding to the kinematic limit x = 1, is determined by 
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the Pauli and Dirac electromagnetic form factors of the nucleon as follows: 



gt\Q 2 ) 



e; = -F 1 (g 2 ) [F.iQ 2 ) + F 2 (Q' 



(16) 



In our numerical analysis we use parametrizations from ref. |[22|| . 



5.1 Resonance contributions 



In order to investigate the role of low-mass nucleon excitations it is useful to introduce 
the ratio 



(n) 
9l 



(Q2) ' J Q 1 dxx^ 1 g 1 (Q 2 



with x = x(W = W ). 



(17) 



In the numerator we sum over all contributions of nucleon resonances and non-resonant 
multi-meson excitations with invariant mass Mx < Wo = 2 GeV. In addition we always 
include the elastic part ([H]). Figure [3| shows that these low mass contributions to g[ n ^ 
are quite sizable, especially for higher moments. For example, at Q 2 = 2 GeV 2 they are 
responsible for about 20% of the first moment , while they account for 75% of gf\ 
With increasing n the role of low-mass excitations becomes evidently more pronounced. 
At the same time, the influence of the low-mass part of the spectrum also increases 
with decreasing Q 2 . Roughly speaking, for Q 2 < 2n GeV 2 low- mass excitations with 
W < 2 GeV account for more than 10% of g[ n \ At large Q 2 continuum contributions 
with W > Wo take over. A similar observation has been made in an analysis of unpolarized 
lepton scattering pfl. 



5.2 Higher twist analysis 

In order to extract the genuine higher twist coefficients f n from the structure function 

(n) 

moments g\ one has to subtract twist-2 contributions and target mass corrections from 
each given moment. Returning to eq.(H) we have: 



g { r\Q 2 )W = gt\Q 2 )-\a n -im-^^^ 

„ AM 2 ^/M 4 \ . , 

= ^w + M' (19) 

In the following we consider the first three moments, n = 1,3,5. We compare results 
obtained from our model for g\ with the twist-2 contributions a n _i/2 from the NLO 



parametrization of ref . ||14|| . We also study the influence of target mass effects. Finally we 

ht- 



discuss different contributions to the higher twist part g[ n \Q 2 ) 
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In Fig.|] we compare the full moments with the leading twist parts, a n _i/2, and 
the higher twist components g± (Q 2 )\u- At small Q 2 one observes significant differences 
between g[ n ^ and a n _i/2. In particular one finds g{ n) {Q 2 )\ ht > 0.1a n _x/2 for Q 2 < 2, 4, 10 
GeV 2 and n — 1, 2, 3, respectively. The region where higher twist becomes important 
depends obviously on the moment n. For fixed Q 2 the difference between g\ L \Q 2 ) and 
a n _i/2 increases with n. This is easily understood since contributions of low- mass nucleon 
excitations are enhanced in higher moments as pointed out in the previous section. Also 
shown in Fig.[| is the size of target mass effects. Since the coefficients d n are not known 
accurately we use d n = which is compatible with present data fl], || and corresponds 
to the Wandzura-Wilczek conjecture p3| . For this choice target mass effects are indeed 



small. As an example, at Q 2 = 2 GeV 2 and n=l,3,5 they amount to less than 10% of the 
higher twist part. To estimate the uncertainty of this result we also use d n obtained from 
eq.@ for g2{x) = 0. In this case target mass effects increase significantly and lead to a 
decrease of <7i (Q 2 )|ht by about 30%. High precision data on the spin structure function 
g 2 , which are of course interesting in their own right, are therefore an important ingredient 
in the QCD analysis of g\ itself. 

Twist-4 contributions to g^ are proportional to 1/Q 2 (up to logarithmic corrections). In 
order to have a closer look at these terms it is instructive to plot the higher-twist moments 
gi l \Q 2 )\ht versus 1/Q 2 , as done in Fig. |^. Their approximately linear behavior indicates 
that twist-4 contributions play indeed a dominant role in g[ n \Q 2 )\ht- Neglecting terms 
of twist-6 and higher gives /| ~ 0.1 at Q 2 = 2 GeV 2 , which agrees quite well with the 
analysis of ref. |5]]. Further estimates can be found in ref.|T^ . 



In the same figure we show the separate contributions to g^iQ 2 )^ from elastic scattering 
and from low-mass excitations with (M + m n ) < W < 2 GeV. Evidently, none of these 
contributions is small, in fact they all are of the same order of magnitude as g[ n \Q 2 )\ht 
itself. 

These observations emphasize the need for high-precision measurements especially in the 
resonance region. Upcoming data from TJNAF are certainly welcome here. Figure 
|5] also points to the crucial role played by the elastic piece ([Ljj). Its proper treatment 
requires accurate information on the nucleon electromagnetic form factors in the range 
1.5 GeV 2 < Q 2 < 10 GeV 2 . 

For the higher moments with n = 3, 5 the kinematic window in which twist-4 contributions 
dominate, that is, where g^iQ 2 )^ behaves linearly with 1/Q 2 , moves successively to 
higher Q 2 . Again the contributions from elastic, resonant and non- resonant scattering all 
turn out to be of similar importance. 
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5.3 Parton-hadron duality 



With decreasing Q 2 the higher twist contributions eventually reach the magnitude of the 
leading twist parts. As a consequence the twist expansion (|J) breaks down. Our model for 
g\ can be used to suggest where this transition takes place: for a given moment g[ n ^ higher 
twist contributions amount to less than 50% of the leading twist ones if Q 2 > n GeV 2 . 



On the other hand, we have learned in section |5.1| that low mass excitations account for 
more than 10% of g[ n) if Q 2 < 2n GeV 2 . 

This indicates a region of n and Q 2 in which perturbative higher twist corrections coexist 
with resonance contributions. The resonance terms are significant, and the transition 
amplitudes involving these resonances introduce powers of 1/Q 2 in just such a way that 
they follow the deep- inelastic, large Q 2 behaviour of g\. 

Such a behavior is known as parton-hadron duality, a notion introduced by Bloom and 
Gilman for the unpolarized structure function F 2 . A QCD explanation of this phenomenon 
has first been offered in ref. [0] and was further elaborated in ref. [^5|. According to our 



results similar arguments apply to polarized lepton-nucleon scattering. 



6 Summary 

i) Contributions from the region of the nucleon resonances are an essential ingredient in 
the "higher-twist" analysis of the spin structure function g\. Their effects are clearly 
visible in g\ even at Q 2 as large as 5 GeV 2 . For example, low mass excitations with 
W < 2 GeV account for more than 50% of the 3rd moment and more than 80% of 
the 5th moment of g\ in the range Q 2 < 3 GeV 2 . 

ii) We have pointed to the importance of the elastic scattering (x — 1) part in a con- 
sistent moment analysis of g±. Without inclusion of this elastic part, an extraction 
of higher-twist terms would not be meaningful. 

iii) We observe a coexistence of resonance contributions and perturbative higher-twist 
corrections in a window, roughly framed by n GeV 2 < Q 2 < 2n GeV 2 , where 
n = 1,3,5, .. . denotes the moment of g\. The understanding of this coexistence 
region in terms of parton-hadron duality is an interesting issue. Precision data from 
TJNAF will help clarifying these questions in the near future. 
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12 3 4 5 

W 2 [GcV 2 ] 

Figure 1: The proton spin structure function g\ at Q 2 = 1.2 GeV 2 as calculated from 
our model in section |j. The data are taken from ref . Jl| . 
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Figure 2: The x-dependence of g\ for Q 2 = 1, 2, 5 and 10 GeV 2 . The full lines show 
results of the model developed in section f|, the dashed lines correspond to the twist- 2 
parametrization of ref.fT4] 
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Figure 3: Contribution to the moments g[ n ^ from low mass excitations with W < W = 2 
GeV, normalized to the full moments. Shown are the 1 st (full), 3 rd (dotted), 5 th (dashed), 
and 7 th (dash-dotted) moments. 
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Figure 4: The moments g[ and their decomposition (|i~8|) for n = 1, 3 and 5. Shown are: 
#i (dashed), the twist-2 part a n _i/2 (dotted), the target mass corrections (dash-dotted), 
and the higher twist piece gi^lht of Eq.(|l8|)|full). 
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Figure 5: Higher twist components gf^lht f° r n — 1,3 and 5 and their decomposition, 
plotted as a function of 1/Q 2 . Shown are gf^lht (full), the elastic part (dotted), contri- 
butions from low- mass excitations with W < W = 2 GeV (dashed), and the difference 
between the high-mass continuum with W > W and the twist- 2 part (dash-dotted). 



17 



